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ABSTRACT 



To shed light onto the circumnuclear environment of 22 GHz (A ~ 1.3 cm) 
H2O maser galaxies, we have analyzed some of their multi-wavelength properties, 
including the far infrared luminosity (FIR), the luminosity of the [OIII]A5007 
emission line, the nuclear X-ray luminosity, and the equivalent width of the neu- 
tral iron Ka emission line (EW (Kq,)). Our statistical analysis includes a total of 
85 sources, most of them harboring an active galactic nucleus (AGN). There are 
strong anti-correlations between EW {Ka) and two "optical thickness parame- 
ters" , i.e. the ratios of the X-ray luminosity versus the presumably more isotrop- 
ically radiated [O III] and far infrared (FIR) luminosities. Based on these anti- 
correlations, a set of quantitative criteria, EW (irQ)>300eV, I/2-iofeey<2 1/foiii] 
and Lfir>600 L2-iofcey can be established for Compton-thick nuclear regions. 
18 H2O maser galaxies belong to this category. There are no obvious correla- 
tions between the EW (K^,), the [OIII] luminosity and the isotropic H2O maser 
luminosity. When comparing samples of Seyfert 2s with and without detected 
H2O maser lines, there seem to exist differences in EW (K^,) and the fraction of 
Compton-thick nuclei. This should be studied further. For AGN masers alone, 
there is no obvious correlation between FIR and H2O maser luminosities. How- 
ever, including masers associated with star forming regions, a linear correlation is 
revealed. Overall, the extragalactic FIR-H2O data agree with the corresponding 
relation for Galactic maser sources, extrapolated by several orders of magnitude 
to higher luminosities. 

Subject headings: Masers - galaxies: active - galaxies: nuclei - galaxies: Compton 
thickness : galaxies - X-rays: galaxies 
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1. Introduction 

Thanks to a large number of dedicated surveys during the last 15 years, the number 
of galaxies known to host 22 GHz (A~1.3cm) H2O masers has increased tenfold and 85 
sources at distances larger than those of the Magellanic Clouds have been reported so 
far to exhibit H2O maser emission (e.g., Braatz & Gugliucci 2008, Darhng et al. 2008, 
Greenhill et al. 2008). Some of them are AG N- related, while others are found in off-nuclear 
star forming regions. Among the masers clearly identified as AGN-related, a large fraction 
(~40%) have been identified as "disk- maser" candidates (Kondratko et al. 2006). Their 
maser spots are associated with central, typically parsec sized molecular accretion disks and 
maser line spectra show "high velocity features" (red-shifted and blue-shifted features), in 
addition to the systemic velocity components. The study of disk-maser sources has become 
a very important subject, permitting mass estimates of supermassive black holes, distance 
estimates of galaxies, and providing a perspective to improve the accuracy of the Hubble 
constant and to constrain the equation of state for the elusive dark energy (e.g. Braatz et 
al. 2009). 

Observations show that among AGN masers, H2O maser spots locate preferentially in 
the nuclear regions of Seyfert 2 or LINER galaxies and most of them are heavily obscured 
(A^H>10^^cm-2; Braatz et al. 1997, Madejski et al. 2006, Zhang et al. 2006, Greenhill et 
al. 2008). The nuclear X-ray source is generally believed to heat the gas to temperatures 
suitable for 22 GHz H2O maser emission (Neufeld et al. 1994), which is supported by a 
relation between maser luminosity, the unabsorbed intrinsic nuclear X-ray luminosity, and 
the mass of the black hole (Kondratko et al. 2006, Su et al. 2008). 

For the obscured nuclear regions H2O masers and X-rays provide, unlike optical data, 
deeply penetrating views. The X-ray absorption along the line-of-sight to the nucleus 
can provide through spectral model fitting important information on the nature of the 
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circumnuclear environment. With modern X-ray telescopes, high quahty X-ray spectra 
have been obtained for more than 30 H2O maser galaxies. Based on an analysis of such 
spectra, the nuclear column densities of maser host galaxies were investigated by Zhang et 
al. (2006) and Greenhill et al. (2008). While these studies have shown that nuclear H2O 
masers are mostly found in environments with high column density, it is still open whether 
H2O masers preferentially arise from Compton-thick (iVH>10^^cm~^) AGN. 

To reduce limitations and uncertainties when modeling X-ray spectra, here we try 
to provide additional constraints to evaluate line-of-sight column densities. The iron 
emission line (~6.4keV) is the most prominent line in the X-ray spectra of AGN. It is 
believed to be produced either by transmission through the absorbing material (Leahy & 
Creighton 1993) or via the process of X-ray scattering/reflection by the cold iron in the 
nuclear accretion disk (e.g., Lightman & White 1988; Fabian et al. 2000) or torus (Ghisellini 
et al. 1994). So it provides useful information to constrain the column density that absorbs 
the continuum and, in some cases, to distinguish between Compton-thick and -thin sources. 
Flat hard X-ray spectra and high equivalent widths (EW (Kq,)) of the iron line were found 
in highly absorbed sources and are generally used to identify Compton-thick nuclei (e.g.. 
Matt 1997). However, currently there are no quantitative criteria on EW (Kq) yet as a 
probe of the gas absorption. 

Unlike the X-ray absorption, [OIII]A5007 and FIR emission have often been used as 
comparatively isotropic indicators of the intrinsic nuclear power (e.g, Mulchaey et al. 1994, 
Alonso-Herrero et al. 1997, Shu et al. 2007). The forbidden [OIII]A5007 line emission 
originates in the narrow line region (NLR), which is assumed to be isotropic for both type 
1 and type 2 systems. Due to possible shielding effects in the torus, Netzer et al. (2006) 
proposed that the [OII]A3727 line is a better tracer than [OIII] for its considerably larger 
emission region. However, this still needs to be checked, since the [O II] line can also 
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arise from high mass star formation. And their results are based on a high-luminosity 
high-redshift AGN sample. For the less luminous, relatively nearby AGN as those of our 
maser sample, the [OIII]A5007 luminosity has been extensively used as an indicator of the 
intrinsic AGN power (e.g., Mulchaey et al. 1994, Alonso-Herrero et al. 1997, Maiolino & 
Rieke 1995, Heckman et al. 2005, Risaliti et al. 1999, Bassani et al. 1999, Panessa et al. 
2006, Lamastra et al. 2009). The FIR emission, produced on linear scales much larger than 
the nuclear torus, should be free of biases caused by the viewing angle and should be highly 
isotropic (e.g., Mulchaey et al. 1994). So even without a detailed knowledge of the X-ray 
spectrum, a comparison of the observed X-ray luminosity with those of the more isotropic 
tracers of nuclear power can help us to determine the gas absorption along the line of sight 
toward the AGN (e.g., Bassani et al. 1999). 

In the following. X-ray, [O III], and FIR data are collected for those H2O maser galaxies 
which are located farther than the Magellanic Clouds and which are known to host 22 GHz 
H2O masers. The data are analyzed to determine gas column densities along the hne of 
sight toward the AGN and to derive properties related to activity in highly obscured and 
therefore particularly elusive Compton-thick nuclear environments. 

2. Data 

Data related to the galaxies with detected H2O maser emission, located at larger 
distances than the Magellanic Clouds, are compiled in Table 1. Among the 85 published 
H2O maser sources, there are 66 AGN masers (63 "megamasers" with isotropic luminosities 
LH2O > 10 L0 and 3 "kilomasers" with LH2O < 10 L0). The megamasers are classified as 
such because of their luminosity. While most of them have not yet been studied in detail, 
all thoroughly investigated megamasers are not separated by more than a few pc from 
the line of sight to the AGN of their parent galaxy. Ten kilomasers are related to star 
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formation regions and another nine kilomasers are still awaiting interferometric observations 
to investigate their nature. All activity types of these H2O maser galaxies are listed in 
Table 1 and their dominant types (following Bennert et al. 2009) are used for the statistics 
(for details, see Table 2). Most of the AGN maser sources are Seyfert 2 galaxies or LINERs. 



- 7- 



Table 1. Physical Parameters of Extragalactic H2O Maser Sources* 



Source 



Type 



Tel. EW (Ka) Fx 



Ref. 



FpiR B.a/Hg F[oiu] 



Ref. 



NGC17 Sy2,LIRG, HII 

NGC23 LINER.LIRG, HI 
ICIO 

NGC 235A Syl 

NGC253 Sy2,SBG,HII 

NGC 262 (Mrk 348) Sy2 



IRAS F01063-803 
NGC449(Mrkl) 
NGC 520 
NGC 598 
NGC 591 
NGC 613 
IC0184 
NGC 1052 
NGC 1068 

NGC 1106 
Mrk 1066 
NGC 1320 
NGC 1386 

IRAS03355+0104 
IC342 

MG J0414+0534 
UGC 3193 

UGC3255 
Mrk 3 



NGC2146 
VII ZW73 



Sy2 

SBG,HII 
HII 

Sy2 

Sy2,HII 
Sy2,HII 
Sy2,LINER 
Sy2,Syl 

Sy2 
Sy2 

Sy2 
Sy2 

Sy2 

Sy2,Hn 
QSOl 

Sy2 
Sy2 



HII 

Sy2 



A 
G 



B 
X 
A 

C 
X 

X 
A 



X 
B 
A 



212 



+68 



230 



+120 

140 



2200 



+700 
-600 



180i«° 
1200 ± 500 
1210+^®" 



2200ll« 
1800t^°° 
7600tS 



610 



+30 



650 
997; 



+182 
182 



+300 
-307 



482 
1270 



20 



400 
462 



23 
496.2 

27 ±5 
20 



590 
650 
650 



AwaOO 
Bas99 



Gua05b 



Ter02 
Cap06 



350 Bas99 



Shu07 

Gre08 
Gua05b 
Bas99 



Bia05 
Cap99 
Bas99 



1.33 
1.09 
3.65 
1.1 
73.56 
0.34 

0.37 
0.32 
2.86 
67.71 
0.28 
2.6 

0.14 
25.01 

0.19 
1.14 
0.29 
0.72 

0.13 
8.29 

0.36 
0.19 
0.57 



12.54 
0.21 



5.89 
12.5 



6.02 



5.86 



4.55 



2.82 
7.00 



37.9 Ho97 
442.85 Ho97 



177 Dah88 

236.51 Dah88 

0.24 Ho97 

178 Whi92 



13.3 Dah88 
6780 Dah88 



1.51 514 Whi92 



5.7 1020 Sto89 



7.69 3.4 Ho97 



6.67 4610 Whi92 



11.1 30.47 Ho97 
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Table 1 — Continued 



Source 



Type 



Tel. EW (Kc.) 



Fx 



Ref. 



Pfir Ha/H^ 



^[OIII] 



Ref. 



NGC 2273 

UGC 3789 

Mrk78 
J0804+3607 
He 2-10 

2MASX J08362280 

Mrk 1210 



NGC 2639 

NGC 2782 

NGC 2824(Mrk394) 

SBS 0927+49 

NGC 2960 

UGC 5101 

NGC 2979 

NGC 2989 

NGCS034 

NGC 3079 

Mrk 34 

NGC3359 
IC 2560 

NGC 3393 



Sy2 



Sy2 

QS02 

SBG 

Sy2 

Sy2,Syl 



Syl.9 

Syl,SBG 

Sy? 

LINER 
LINER 

Syl.5,LINE R.LIRG 

Sy2 

HII 

SBG,HII 
Sy2,LINER 

Sy2 
HII 
Sy2 

Sy2 



NGC 3556 HII 
ARP299(NGC 3690) 



X 
A 



C 
X 
B 
A 
A 
C 



X 
B 



X 
C 
X 
A 

X,C 
B 



2200+**°° 



1040 



69 



+440 
-460 



~188 
130 



108 



+50 
-65 



8201^60 
14901111^0 
990 



410 



+270 
-240 



1480 ± 500 
2400t?§°° 



2320lifo 
2770+490 
1400 ± 800 
3500 ± 2000 



^22+262 
*^^-288 



636 



+236 



12 

125 



9.56 



+0.63 
-0.64 



840 
970 
930 
160 
25.3 
~30 



33 
37±8 



38.i 
38.4 



1-1.8 
-5.1 

+21.1 



□+6 



40 



43.7 



Gua05b 
Ter02 



Jur05 

Zha09 
Gua02 
Ohn04 
AwaOO 
Ter02 
Zha06 



Ima03 



Cap06 
lyoOl 



Til08 
Mad06 
Gu05a 
Bas99 

Bal04 
Del02 



0.7 

0.26 
0.15 

2.40 

0.17 
0.36 



NGC 3735 



Sy2 



0.43 
0.91 
0.12 
0.29 
0.24 
1.06 
0.2 

96.49 
5.5 

0.43 
0.70 
4.29 

0.33 

7.29 
9.43 

1.02 



6.92 



6.50 



5.20 



4.16 
6.67 



25.0 
25.0 

10.5 



4.12 

7.14 
5.88 

6.31 



164 Whi92 



242.76 Dah88 



580 Ter91 



4 Ris99 
62.19 Ho97 



191 Kim95 

1615.17 Ho97 

92 Ho97 

204.19 Dah88 

> 40 Ris99 

316 Dia88 

2.24 Ho97 

35.6 Ho97 

33 Ho97 



Table 1 — Continued 



Source 


Type 


Tel. 


EW {Kc) 


Fx 


Ref. 


Ffir 


H„/H^ 


-P[OIII] 


Ref. 


Antennae 


SBG 


C 
















NGC 4051 


Syl.5 


X 


240 ± 40 


627 


Cap06 


1.32 


3.33 


59.99 


Ho97 


NGC4151 


Syl.5 


X 


300 ± 30 


4510 


Cap06 


1.11 


3.45 


1695.39 


Ho97 






A 


101±5 


~20000 


WeaOl 










NGC4214 


SBG 


C 




243 


Har04 


1.81 








NGC 4258 


Syl.9,LINER 


X 


27 ±20 


837 


Cap06 


5.08 


9.12 


262 


Hec80 






A 


250±100 


300 


Bas99 










NGC 4293 


LINER 










0.58 


7.69 


5.95 


Ho97 


NGC 4388 


Sy2 


X 


440 ± 90 


762 


Cap06 


1.44 


5.50 


374±50 


Bas99 






A 


709+243 


1200 


Bas99 










NGC 4527 


LINER,HII 










3.80 








ESO 269-G012 


Sy2 










0.23 








NGC 4922 


Sy2,LINER 










0.58 


7.14 


33.03 


Kim95 


NGC 4945 


Sy2 


C 


1300 


500 


Don03 


41.38 




> 40 


Ris99 






B 


~1300 


540 


GuaOO 














A 


850 ± 160 


350 


Bas99 










NGC 5194 


Sy2,HII 


X 


986 ±210 


48 


Cap06 


6.62 


8.33 


228 


Ho97 






A 


^"-360 


91.9 


Ter02 










NGC 5253 


SBG,HII 


C 




29.9 


Jur05 


3.16 








Mrk266(NGC 5256) 


Sy2,LIRG,SBG 


B 


575 


56 


RisOO 


0.84 


5.92 


44.33 


Dah88 


NGC 5347 


Sy2 


C 


1300 ± 500 


22 


Lev06 


0.27 




114 


TraOl 


NGC 5495 


Sy2,HII? 










0.27 








Circinus 


Sy2 


C 




1400 


SmiOl 


26.11 


19.1 


6970 


Bas99 


NGC 5506(Mrk 1376) 


Syl.9 


X 




5800 


Bia03 


1.06 


7.20 


333 


LumOl 






A 


150 ± 30 


8380 


Bas99 










NGC 5643 


Sy2 


X 


500 


84 


Gua04 


2.59 


6.40 


662 


Whi92 






A 


18001^0 


130 


Bas99 










NGC 5728 


Sy2,HII 


C 




133 


Shu07 


0.96 


5.96 


761 


Sto95 






C 


1130 




Zha06 










UGC09618NED02 


LINER,HII 










0.92 








NGC 5793 


Sy2 










0.6 








NGC 6240 


Sy2,LINER 


C 


2400lf°° 


170 


Pta03 


2.18 


17.2 


135±20 


Kim95 



Table 1 — Continued 



Source 


Type 


Tel. 


EW (K^) 


Fx 


Ref. 


Ffir 




-'^[OIII] 


Ref. 






A 


15801^0 


190 


Bas99 










NGC6264 


Sy2 


















NGC 6300 


Sy2 


X 


148 ± 18 


2160 ± 100 


Mat04b 


2.34 




320 


LumOl 


NGC6323 


Sy2 


















ESQ 103-G035 


Sy2,Syl 


A 


1 7Q+60 


907 


Tur97 


0.37 


6.31 


112 


Pol96 


IRAS F19370-013 


Sy2,HII 










0.26 








3C403 


FRII 


C 


244±20 




Kra05 


0.1 








NGC 6926 


Sy2,HII 










0.65 


15.6 


21.98 


Kim95 


AM 2158-380NED02 Sy2,RG 


















TXS 2226-184 


LINER 


















NGC 7479 


Sy2,LINER 










0.83 


10.0 


37.9 


Kim95 


IC 1481 


LINER 










0.12 









Note. — 85 published extragalactic H2O maser sources with available physical parameters are listed (10 masers arise in star 

forming regions marked by italics and 29 out of 66 AGN-masers are potential disk-masers, their source names are presented in 
boldface). Recently reported masers, not being part of the 78 sources listed by Bennert et al. (2009), are one type I quasar 
from Impellizzeri et al. (2008), two sources (NGC 17 and NGC 1320) from Greenhill et al. (2008), and four new masers related 
with star formation (He 2-10, Antennae, NGC 4214, NGC 5253) from Darling et al. (2008). 

^'Extragalactic H2O maser host galaxies; 

Coi«mn2:rpypg of nuclear activity SBC; StarBurst Galaxy; Syl, Syl.5, Syl.9, Sy2: Seyfert types; LINER: Low-Ionization 
Nuclear Emission Line Region; LIRG: Luminous-Infrared Galaxy; FRII: Fanarov- Riley Type II radio galaxy; NLRG: Narrow- 
Line Radio Galaxy; RG: Radio Galaxy; HII: classified as a H II region; QSOl and QS02: type 1 and 2 Quasars. References: 
Zhang et al. (2006); Kondratko et al. (2006) and NED; 

CoiMmn3:x-ray telescope-A: ASCA; B: BeppoSAX; C: Chandra; X: XMM-Newton; 

Coi«mn3 4&B:rpjjg (j^^^ ^.^^ [j^^^ (^y) ^.^^ 2-10 keV observed X-ray flux (in units of 10"" ergs'^cm-^); 

6: References for Col. 4&5 — AwaOO: Awaki et al. 2000; Bal04: Ballo et al. 2004; Bas99: Bassani et al. 1999; Bec04: 
Beckmann et al. 2004; Bia03: Bianchi et al. 2003; Bia05: Bianchi et al 2005; Cap99: Cappi et al. 1999; Cap06: Cappi et al. 
2006; Del02: Delia Ceca et al. 2002; Dia88: Diaz et al. 1988; Don03: Done et al. 2003; GuaOO: Guainazzi et al. 2000a; Gua02: 
Guainaazi et al. 2002; Gua04: Guainaazi et al. 2004; Gua05a: Guainaazi et al. 2005a; Gua05b: Guainazzi et al. 2005b; Ima03: 
Imanishi et al. 2003; Iwa02: Iwasawa et al. 2002; lyoOl: lyomoto et al. 2001; Jen04: Jenkins et al. 2004; Kra05: Kraft et 
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al. 2005; Lev06: Levenson et al. 2006; Mad06: Madejski et al. 2006; MatOl: Matt et al. 2001; Mat04a: Matt et al. 2004a; 

Mat04b: Matsumoto et al. 2004b; Ohn04: Ohno et al. 2004; Pta03: Ptak et al. 2003; RisOO: Risaliti et al. 2000; Smi96: Smith 
& Done 1996; SmiOl: Smith & Wilson 2001; Ter02: Terashima et al. 2002; Tur97: Turner et al. 1997; WeaOl: Weaver et al. 
2001; Zha06: Zhang ct al. 2006; Zha09: Zhang et al. 2009; 

Cohirnn 7: pjj^ jj^ ^j^j^g gj-gs"! Cm-^; 

Column 8:rpjjg Balmer increment line intensity ratio Ha/H^j; 

Coiumn 9: Extinction corrected [OIII]A5007 line flux in units of 10"" erg s'^ cm.-'^; 

Coiumn 10: j^gfgrences for Col. 8&9 — Bas99: Bassani et al. 1999; Dah88: Dahari & De Robertis 1988; Hec80: Heckman et 

al. 1980; Ho97: Ho ct al. 1997; Kim95: Kim ct al. 1995; LumOl: Lumsden et al. 2001; Pol96: PoUctta ct al. 1996; Ris99: 
Risaliti ct al. 1999; Shu07: Shu ct al. 2007; Sto89; Storchi-Bcrgmann & Pastoriza 1989; Sto95: Storchi-Bergmann et al. 1995; 
Ter91: Terlevich et al. 1991; TraOl: Tran 2001; Whi92: Whittle 1992; 
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Besides the type of nuclear activity, we collected multi-wavelength data and parameters 
for all of the 85 extragalactic H2O maser sources, including the EW (Ka) of the iron 
emission line, the observed X-ray (2-10 keV) flux Fx, the FIR flux Fpui, the observed 
Ra/Hp line intensity ratio and the [OIII]A5007 flux density, Fjom]. These data are also 
presented in Table 1. 

The X-ray data of our sample are collected from the literature, based on observations 
from ASCA, BeppoSAX, Chandra, and the XMM-Newton satellite. For some sources more 
than one value were reported in the literature for parameters such as the EW (K^) of 
the iron line and the observed X-ray flux, either due to intrinsic variability of the source 
or due to a different modeling of the spectra. For comparison, all results available from 
the literature are listed in Table 1. For our statistical analyses, observational results were 
taken from XMM-Newton and Chandra whenever possible. Otherwise the most recent 
measurements were used. 

For the [OIII]A5007 emission line, the fluxes are mainly compiled from Dahari & De 
Robertis (1988), Ho et al. (1997), Kim et al. (1995) and Whittle (1992). Again the most 
recent result was adopted when two or more observations were found. Since the NLR may be 
obscured by material in the host galaxy, the observed [OIII]A5007 flux should be corrected 
to determine its intrinsic flux. The narrow-line Balmer decrement was used to estimate this 
extinction. Assuming an intrinsic Balmer decrement {Hoc/Hp)^ — 3, the intrinsic [OIII] 
line fluxes were derived from the formula F[oiii],cor = F^om],ohs[{Ha/ Hp)obs/ {Ha/ Hp)^]^-^'^ 
(Bassani et al. 1999). The [OIII]A5007 fluxes were obtained for 46 maser sources and their 
luminosities were calculated from their extinction-corrected fluxes. 

The IRAS (Infrared Astronomical Satellite) point source catalogue was used to obtain 
the infrared fluxes (12, 25, 60, 100//m) for our H2O maser host galaxies. Following the 
method of Wouterloot & Walmsley (1986), the infrared flux (6 < A < 400//m) was derived 
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by extrapolating flux densities beyond 12 and 100//m and assuming a grain emissivity 
proportional to frequency v. In this way, infrared luminosities could be determined for 
76 H2O maser galaxies. Throughout the paper, the luminosity distance was derived using 
Calculators I provided by the NASA Extragalactic Database (NED), assuming = 0.270, 
= 0.730, and //o=70kms-^ Mpc"^ (e.g., Spergel et al. 2007). 

3. Results 

3.1. Isotropic indicators of the intrinsic nuclear power 

For our H2O maser galaxies with available data, the luminosity distributions of both 
isotropic indicators, the [OIII] line and the FIR luminosities, are presented in Fig. 1. The 
left panel shows the histogram of the infrared luminosities for all 76 H2O maser galaxies 
with available data {logLpm, hereafter luminosity in logarithmic scale and in units of 
ergs~^). Log LpiR ranges from 41.8 to 45.5, with a mean value of 44.1±0.1 (the error 
denotes the standard deviation of the mean). For the sub-sample of 10 masers associated 
with massive star formation, the mean value of <log Lfir>—^?). 2)^^.2 is slightly fainter 
than that of the AGN masers. However, the difference between both distributions (not 
shown here) is not significant. The histogram shown in the right panel presents the number 
of the H2O maser sources as a function of the [OIII] hue luminosity. It gives the range 
of <Zog'L[oiii]> from 39 to 43.2 and a mean value of 40.9±0.3, comparable to results of 
Seyfert 2 galaxies not associated with detected H2O maser (Mulchaey et al. 1994). For 
comparison, the latter sample is also presented in Fig. 1 (dashed lines). 

For the two tracers Log Lpm and Lo^fL^oiii], the ratio was calculated and the 
distribution of the ratio is shown as a histogram in Fig. 2 (left panel). No significant 
difference can be found for the distribution of the ratio between our H2O maser galaxies 
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Fig. 1. — Distributions of the FIR luminosity {logLpiR, 76 sources, left panel) and the 
[OIII]A5007 line luminosity (log [Ojjj], 46 sources, right panel) for H2O maser sources (solid 
lines), in units of ergs~^. For comparison, the distributions are also given for a Seyfert 2 
sample without known H2O maser emission (dashed lines, 38 sources from Mulchaey et al. 
1994). All numbers are plotted on a percent scale (%). 
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Fig. 2. — Left panel: The distribution (on a percent scale) of the ratio between the FIR and 
the [O III] line luminosity, corrected for extinction, of H2O maser galaxies on a logarithmic 
scale. For comparison, the distribution is also given for a Seyfert 2 sample without known 
H2O maser emission (dashed lines, 38 sources from Mulchaey et al. 1994). Right panel: the 
FIR versus extinction-corrected [O III] luminosity. The straight line shows a linear fit to the 
unweighted data. 
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and the Seyfert 2 galaxy sample not containing galaxies with known masers (Mulchaey 
et al. 1994). This is in agreement with the fact that most H2O maser sources have been 
found in Seyfert 2 systems. The good agreement between Seyfert 1 and Seyfert 2 galaxies 
(Mulchaey et al. 1994) with respect to this parameter, in spite of different viewing angles, 
lends further support to an approximately isotropic emission of both tracers. In order 
to further compare those two isotropic tracers, the FIR luminosity is plotted against the 
extinction-corrected [O III] luminosity of our maser sample in the right panel of Fig. 2. A 
correlation is found, logLi?7i^=(30.86±2.30)+ (0.32±0.06)logL[om] with a Spearman's rank 
correlation coefficient r=0.66 and a chance probability < 5xl0~^. Assuming the corrected 
[OIII] line luminosity to be a good isotropic tracer (Bassani et al. 1999), the strong 
correlation between FIR luminosity and [O III] line luminosity may suggest that the FIR 
luminosity is to some extent also an indicator of intrinsic nuclear activity, even though the 
FIR flux may be contaminated by a spatially extended starburst component. 

3.2. Optical thickness parameters and the 6.4 keV iron line 

As mentioned above, nuclear absorbing columns can be obtained by analyzing X-ray 
spectra. These analyses are model dependent. Comparing the observed 2-10 keV X-ray 
emission from the nuclear region, absorbed by the obscuring material along the line-of-sight, 
with the intrinsic nuclear power provides another method for cvahiating the absorbing 
column density. The ratio of the observed 2-lOkeV luminosity and the luminosity of the 
nuclear isotropic indicators was assumed to represent the optical thickness parameter, 
which allows us to diagnose the gas absorption (Bassani et al. 1999). Here the Lx/L[oni] 
and hx/^FiR ratios are used as optical thickness parameters. In addition, high EW {Ka) 
values of the iron emission line are considered as a qualitative feature indicating a heavily 
obscured nucleus (e.g., Maiohno & Risaliti 2007). Combining the EW {Ka) and the optical 
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Fig. 3. — EW {Ka) of the iron line versus the thickness parameters L2-ioitey/L[oiii] (Tl, 
left panel) and ^2-wkeV I^fir (T2, right panel). Square and pentacle symbols represent 
Compton-thick and -thin sources as identified by modeling their X-ray spectra (e.g., Zhang 
et al. 2006). The bold squares show those Compton-thick sources, which have high EW 
(Kq) and low T values. These are located in the upper-left regions of both panels. Three 
ea;ce;?iiona/ sources NGC2639, MrkS, and NGC4945 are marked. For NGC4945, the lower 
limit of its [OIII] luminosity is taken. Details are given in Sect. 3.2. 
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thickness parameters, we probe the circumnuclear environment of the maser sources. 

For the maser galaxies with available data (EW {Ka), Lx, L[oiii] and Lpjn, in total 
31 AGN maser sources), the EW (Kq) of the iron lines is plotted against the optical 
thickness parameters in Fig. 3. The prominent feature of the figure is the existence of an 
ant i- correlation, which is similar to that obtained for the Seyfert 2 sample of Bassani et al. 
(1999), which is predominantly containing targets without known maser lines. In the left 
panel (EW (Ka) v.s. Lx/L[oiii]), a least-squares fit shows log(EW (i^a))=(2.74±0.08)+ 
(-0.31±0.08)logLx/L[oiii], with Spearman's rank correlation coefficient R=-0.57 and a 
chance probability of P ~0.001. The possibly Compton-thick sources (shown by squares in 
Fig. 3), determined by modeling their X-ray spectra, cluster in the upper left region, the 
region with high EW (K^) and low Lx/L[oiii] values. This area is marked by dashed lines. 
The approximate boundaries between Compton-thick and -thin sources, the latter shown 
by pentacles, are log (Lx/L[om])~0.25 (i.e., iv2-iofeey~2 L[om]) and logEW (X„)~2.5 
(i.e, EW~300eV). For comparison with Seyfert 2 galaxies without detected H2O maser 
emission, the Seyfert 2 sample of Bassani et al. (1999) is used, excluding those objects with 
maser emission. The trend is the same for this sample (here not shown). Compton-thick 
sources fall again into the upper left region with high EW (Kq,) and low Lx/Ljom] values, 
while Compton-thin sources are located in the lower right. The right panel of Fig. 3 shows 
a similarly clear anti-correlation between EW (K^,) and Lx/I-'fir- Linear fitting results in 
log(EW {K^))^{1.71 ± 0.26) + (-0.36± 0.08)\ogLx/h fir, with R=-0.61 and P ~ 2 xlQ-^ 
Most of those sources with high EW (Kq,) and low Lx/hpiR value (upper- left region, in 
square symbols) are Compton-thick sources, while those sources with low EW (K^) and 
high Lx/l-'FiR value (in pentacles) are Compton-thin as determined from X-ray spectral 
fitting. The approximate boundaries between Compton-thick and -thin environments are 
in this case logLx/Li;'7R~-2.75 (i.e., LFiR~600L2_iojtey) and logEW {Ka)^2.5, i.e. again 
EW-300eV. 
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We like to emphasize that these boundaries are not arbitrary. Among the 20 Compton- 
thick candidate sources from our sample, classified by conventional X-ray spectroscopy, we 
find 19 in each of our limited "Compton-thick" regions, related either to Lx/L[oiii] or to 
hx/l^FiR- 18 of these sources are identical (see also Sect. 4). 

The results obtained so far could be affected by systematic errors in the measurements. 
Assuming for the [O III] values uncertainties of ~20% (Dahari & De Robertis 1988) results 
in Lx/L[oiii] errors of order 0.08 dex. Obviously, this does not affect the robust fundamental 
trend in our diagnostic Fig. 3. Compton-thick sources are still placed in the upper left and 
Compton-thin sources in the lower right. The anti-correlation between the EW {K^) and 
the optical thickness parameters are readily explained. With an increase of the absorbing 
column density, the X-ray luminosity will decrease so that Lx/L[oiii] and \jx/^fir are 
reduced with respect to Lpiii] and hpiR- On the other hand, the EW {Ka) values will 
increase, since these are measured against a reduced 6.4 keV continuum level. 

With Compton-thick galaxies being located in the upper left part of the panels in Fig. 
3, we find for our H2O maser galaxies three criteria hinting at a Compton-thick nuclear 
environment: EW (irQ,)>300eV, L2_iofeey<2 -C/ioiii] ^-nd I/fir >600 Lo 

-lokev- These are 

independent of the detailed shape of the X-ray spectrum. 

3.3. Indicators for H2O maser emission 

The search for new extragalactic H2O masers is ongoing and important with respect to 
several key aspects of modern astrophysics (see, e.g.. Sect. 1 and Braatz et al. 2009). Here 
we therefore investigate possible indicators of H2O maser emission. We analyze relationships 
between H2O maser luminosity and the iron line EW (K^), FIR, and [O III] luminosity. 

For 33 maser sources with available EW (Ka) of the iron fine in Table 1, we obtain 
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Table 2: Activity types of H2O maser host galaxies (79 sources with available type) 



Type") 


SF-masers'') 


Kilomasers of unknown origin'^' 


AGN-masers'^) 


the whole sample 


Seyfert 2 





4 


45 


49 


LINER 





2 


8 


10 


Inter Sy. 





1 


2 


3 


Seyfert 1 








2 


2 


SBG 


5 


1 





6 


HII 


4 


1 


1 


6 


FRII 








1 


1 


Quasar 








2 


2 



a) Types from NED: Inter Sy., intermediate Seyfert types; LINER, low-ionization nuclear emission line region; SBG, StarBurst 
Galaxy; FRII, Fanaroff- Riley type II radio galaxy; HII, classified as HII region galaxy; b) SF-masers, H2O masers associated 
with off-nuclear star formation regions; c) kilomasers without known type; d) AGN-masers, H2O masers associated with AGN. 
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Fig. 4. — Distributions (on a percent scale) of the iron line EW (Ka) (logarithmic scale, in 
units of eV) for Seyfert 2s with detected H2O maser emission (19 sources, solid lines) and 
non-masing Seyfert 2s (34 sources from Bassani et al. 1999, dashed lines). 
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an EW (Kq) mean value of 945 ± 135 eV. X-ray observations show that strong iron hne 
emission is common in the spectra of AGN hosting H2O masers. It is interesting to check 
if the 6.4 keV hne can be used as a criterion to search for AGN masers. Have Seyfert 2s 
with detected maser emission higher EW (K^) values than non-maser Seyfert 2s? Our 
statistical results show that the iron line EW (Ka) of masing Seyfert 2 galaxies (mean value 
1063 ± 169 eV and median value ~800 eV for our 19 mascr sources) is higher than that of 
the non-maser Seyfert 2 sample (mean value 375 ± 60 eV and median value ~200eV for 
34 sources from Bassani et al. 1999 ). Figure 4 shows the distributions of both samples. 
While the difference seems to be obvious at first sight, we should cautiously avoid a 
definite conclusion due to the large scatter, the still too small number of sources, and 
the incompleteness of the studied samples. Potential differences in sensitivity have also 
to be addressed. For our maser host Seyfert 2s, the EW (Ka) values were taken almost 
exclusively (except four sources) from XMM-Newton or Chandra observations. The results 
for Seyfert 2s without known masers (Bassani's sample) come mostly from ASCA data of 
lower sensitivity, which might lead to an increase in the real average (only sources with a 
strong iron line could be detected). This strengthens our result and amplifies the difference 
between Seyfert 2s with and without detected 22 GHz H2O maser. Nevertheless, we consider 
our result as tentative. 

In order to investigate possible correlations between H2O maser and iron emission 
lines, the EW (Ka) of the Fe line was plotted against the isotropic H2O luminosity in Fig. 5 
(upper panel). As already mentioned, H2O maser emission can be produced by coUisional 
pumping in a dense molecular layer, which is heated by irradiated X-rays from the nucleus 
(Neufeld et al. 1994). Strong Fe emission is believed to be produced via X-ray reflection 
by the cold iron in the circumnuclear region (e.g., Fabian et al. 2000). In those cases where 
H2O maser and iron line emission are detected, the nuclear X-ray emission plays a key role. 
Correlations between maser and iron line emission are therefore expected. However, our 
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Fig. 5. — Upper panel: the EW {Ka) of the iron hne (logarithmic scale, in units of eV) versus 
isotropic luminosity of the H2O maser emission (logarithmic scale, in Lq). Center: [O III] line 
luminosity against Lh^q] Bottom: FIR luminosity against L-h^q, empty circles show Galactic 
H2O masers from Jaffe et al. (1981) and the line marks the correlation of Lh2o/-^[fir]~10^^. 
Squares and triangles represent AGN-masers and masers in star formation regions respec- 
tively. Disk-masers, as a subsample of AGN-masers, for which maser emission is anticipated 
to be particularly well connected with indicators of the intrinsic nuclear power, are marked 
by crosses over squares. 
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results show no apparent trend between the iron hne EW (K^) and the isotropic H2O maser 
luminosity, neither for the entire sample nor the two subsamples, AGN- and star-forming 
masers (squares and triangles, respectively, in Fig. 5). In view of alternative H2O excitation 
mechanisms (see, e.g., Lo 2005), the subsample of possible disk-masers (with detected high 
velocity maser features) was analyzed separately (see the crosses in Fig. 5). However, even 
in this case no significant correlation can be found. 

The [O III] and FIR luminosities were also plotted against H2O maser luminosity in 
Fig. 5 (central and bottom panel respectively). There is no significant correlation between 
L[oiii] and L^^o, although maser sources related to star formation seem to have lower 
[O III] luminosities than AGN masers. For FIR versus H2O maser luminosities, there 
appears to be a correlation, similar to that previously found by Henkel et al. (2005), 
Castangia et al. (2008), Bennert et al. (2009) and Surcis et al. (2009). The relation 
was first found for galactic star forming regions by Jaffe et al. 1981, i.e., luminous H2O 
masers form in star formation regions with high FIR luminosity. For comparison, values 
of Galactic H2O masers are also plotted in Fig. 5 (empty circles) and the line shows 
the correlation of LH2o/-f'[FiR]~10~^ from Jaffe et al. (1981). Apparently, there exists a 
correlation between FIR and H2O maser luminosity over many orders of magnitude. When 
considering AGN-masers only (squares in Fig. 5, lowest panel), the strongest masers appear 
to be overluminous with respect to the Lfir-Lh^q correlation. This is likely caused by the 
different properties of AGN versus star-forming masers. 

4. Discussion 

Compton-thick nuclei are known to contribute a significant fraction of the hard X-ray 
background. Their density as a function of redshift is also a relevant parameter when 
studying the evolution of the universe. Thus our newly defined criteria identifying such 
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nuclei may be helpful when trying to reach this goal. 

Based on our new approach introduced in Sect. 3.2, 18 H2O maser sources in the 
upper left part of the panels in Fig. 3 (see the sohd squares) can be considered to 
be Compton-thick. These are NGC591, NGC1068, Mrkl066, NGC1386, NGC2273, 
UGC5101, NGC2782, NGC3079, IC2560, NGC3393, Arp299, NGC5194, Mrk266, 
NGC 5347, Circinus, NGC 5643, NGC 5728 and NGC 6240. Comparing this with Zhang et 
al. (2006), there are five new sources, NGC 591, UGC5101, NGC 3072, IC2560, Mrk266. 
About 60% of the 31 AGN masers turn out to be Compton-thick. This is consistent with 
the result found by Greenhill et al. (2008). For maser sources associated with Seyfert 2 
nuclei, 60% (15/25) are Compton-thick, which is, however, not significantly, higher than 
that of Seyfert 2 objects without detected maser emission (~45%, 9/20 from Risaliti et al. 
1999; 19/42 from Bassani et al. 1999). 

Our sample also contains a few exceptional sources, which show the limits of our 
selection criteria. For NGC 4945, its low [OIII] line flux places the source outside our 
Compton-thick border (EW {Ka)-Lx/L[om]: see left panel in Fig. 3). However, the [OIII] 
flux from Risaliti et al. (1999) only gives a very stringent lower limit to the actual intrinsic 
[O III] emission. The lack of a reliable [O III] flux is thought to be due to high absorption 
in the edge-on galactic disk, instead of its intrinsic weakness. This is supported by its 
hard X-ray spectrum, which indicates that NGC 4945 hosts one of the brightest AGN in 
the hard X-ray range (>100keV) and can therefore be considered to contain a 'bona-flde' 
Compton-thick nucleus (e.g., Guainazzi et al. 2000b). Similar to NGC 4945, Mrk3 is also 
considered to be a 'Bona-flde' Compton-thick Seyfert 2 from its large brightness in the hard 
X-ray range (e.g., Cappi et al. 1999). However, the source is found outside our hmited EW 
{Ka)-L'x/^FiR region (right panel in Fig. 3) for its relatively low value of the infrared flux. 
More constraints are desirable to probe its circumnuclear environment. The disk- maser 
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galaxy NGC2639 is located inside the required EW {Ka)-Lx/^FiR region and outside 
the EW {Ka)-Lx/^[oiii] area. However, large uncertainties in the X-ray results (ASCA 
observations) have to be noted. NGC 2639 is a weak X-ray source with an ASCA count rate 
<0. 01 counts s~^. The ASCA data were analyzed by Wilson et al. (1998) and Terashima et 
al. (2002) and the source was considered to be Compton-thin (A^h~5x10 ^^cm~^). Low 
photon statistics lead to uncertainties in the fitting models and Chandra observation are 
therefore needed to investigate its highly obscured nucleus. 

We conclude with some cautionary notes. First, uncertainties are involved when using 
the [OIII]A5007 line and FIR luminosity as nuclear isotropic indicators. Although assumed 
to be isotropic, the [O III] line luminosity might depend on the geometry of the system, 
for example on the opening angle of the torus and the inclination of the large scale disk. 
Shielding effects may affect the ionizing radiation seen by the NLR. FIR emission of H2O 
megamaser galaxies may mainly arise from the AGN environment, but there is possible 
contamination from a starburst component and it is unclear how much it contributes. 
Second, the EW (K^) of the iron line can be affected by other factors, such as the geometry 
of the accretion disk and the inclination angle at which the reflecting surface is viewed (e.g., 
Fabian et al. 2000 , Bianchi et al. 2005). A high EW (Kq) of the iron hue can also appear 
if the radiation is anisotropic or if there is a time lag between a drop in the continuum and 
the line emission (Bassani et al. 1999). With future advanced X-ray telescopes, sensitive 
observations of more H2O maser host galaxies at higher energies (above 10 keV) will provide 
important complementary information, further constraining nuclear column density. 



5. Summciry 



In this paper, multi- wavelength data from the complete sample of galaxies (£)>100kpc) 
so far reported to host 22 GHz H2O masers are analyzed, including the equivalent width 
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(EW) of the iron line, the [OIII]A5007 hne, and the X-ray and FIR emission. The latter 
two are considered to be the isotropic tracers of intrinsic nuclear power. The observed 
nuclear X-ray luminosity, compared with the luminosities of these two isotropic tracers (our 
optical thickness parameters), can be used as a measure of the circumnuclear absorption. 
The EW {Ka) and the optical thickness parameters (Lx/I/[oiii], Lx/Lfir) are combined 
here to probe the obscuration of maser host AGN. The main results are summarized below: 

(1) Our statistical analysis shows obvious anti-correlations between the EW {Ka) of the 
Fe emission line and the two optical thickness parameters. Without requiring a full X-ray 
spectrum, Compton-thick nuclear environments can be identified with these parameters 
and are found to be characterized approximately by EW (ira)>300eV, L2_iofeey<2 -L[oiii] 
and LFiR>600L2-iofcev; 

(2) 18 H2O maser sources matching these criteria are identified to be Compton-thick. A 
comparison with Zhang et al. (2006) shows, that among these there are five newly identified 
H2O maser galaxies which are Compton-thick, i.e., NGC591, UGC5101, NGC3072, IC2560 
and Mrk 266. Masers associated with Seyfert 2 nuclei may be more hkely Compton-thick 
(60%) than Seyfert 2s without detected maser emission (~45%). However, this still has to 
be confirmed with larger samples; 

(3) In an attempt to guide future H2O maser surveys, new ways to find extragalactic 
H2O sources are also explored. H2O maser sources may show larger EW (K^) values than 
non-maser Seyfert 2s, which, however, also needs further support. No significant correlations 
have been found between EW (K^), Ljom] and L-^^o- There appears a linear correlation 
between LpiR and I/H2O, which is consistent with the correlation found for Galactic H2O 
masers. However, the strongest H2O masers appear overluminous with respect to their 
LpiR- This may be related to their different origin when compared with masers associated 
with sites of massive star formation well outside the nuclear region of their parent galaxy. 
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